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Abstract 

We compute the quark and antiquark transversity distributions in the nucleon at a 
low normalization point (//~600MeV) in the large-iV c limit, where the nucleon can 
be described as a soliton of an effective chiral theory (chiral quark-soliton model). 
The flavor-nonsinglet distributions, 5u(x)—5d(x) and 5u(x)— Sd(x), appear in leading 
order of the l/iV c -expansion, while the flavor-singlet distributions, Su(x) + 5d{x) 
and 5u(x) + Sd(x), are non-zero only in next -to-leading order. The transversity 
quark and antiquark distributions are found to be significantly different from the 
longitudinally polarized distributions Au(x)±Ad(x) and Au(x)±Ad(x) , respectively, 
in contrast to the prediction of the naive non-relativistic quark model. We show that 
this affects the predictions for the spin asymmetries in Drell-Yan pair production 
in transversely polarized pp and pp collisions. 
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1 Introduction 



A central property of QCD is the possibility to factorize the cross sections for some hard 
scattering processes involving hadrons into the cross section for the partonic subprocess 
involving quarks and gluons, calculable within perturbative QCD, and certain functions 
describing the transition from hadrons to quark and gluon degrees of freedom. Of the 
latter only the scale dependence can be predicted from perturbative QCD; the value 
of the functions themselves at a low scale can only be inferred from experiment, or be 
calculated using non-perturbative methods. In inclusive hard scattering off the nucleon 
[deep-inelastic scattering (DIS)], semi-inclusive particle production, and Drell-Yan pair 
production, the relevant characteristics of the nucleon are the so-called quark and an- 
tiquark distributions in the nucleon - functions describing the emission in the collinear 
direction and subsequent absorption by the nucleon of a quark/ ant iquark carrying a frac- 
tion x of the nucleon momentum. At twist-2 level, helicity conservation allows for three 
different types of such functions Q. The usual unpolarized and longitudinally polarized 
distributions, q(x), q(x) and Aq(x), Aq(x), correspond to probabilities for emitting and 
absorbing a quark or antiquark of same helicity, and can be interpreted as the sum and 
difference of probabilities to find a quark/antiquark with longitudinal polarization parallel 
or antiparallel to that of the nucleon. The third kind of functions, Sq(x), Sq(x), describe 
the emission and absorption of quarks/antiquarks with different helicities [[|. They can 
be interpreted as the probability to find a transversely polarized quark/antiquark in a 
transversely polarized nucleon and are therefore referred to as transversity distributions. 

The crucial difference between the transversity and the "usual" distributions lies in 
their chirality properties. The unpolarized and longitudinally polarized distributions pa- 
rameterize the chirally even parts of the quark density matrix in the nucleon (Dirac struc- 
tures 7 + and 7 + 7s, respectively). They can therefore be measured in DIS at leading-twist 
level, where the chirality of the quarks is preserved by the hard scattering process. These 
distributions are by now well-known, with most of the information coming from QCD fits 
to DIS data |, |5|, In contrast, the transversity distributions describe the chirally 
odd part of the quark density matrix (Dirac structure a +± -y^,). Consequently, they can 
be measured only in hard processes where they enter together with other chirally-odd 
objects, chirally-odd distribution or fragmentation functions. A variety of such processes 
- both hadron-hadron and lepton-hadron induced - are currently being considered; see 
Ref. for a review. In hadron-hadron collisions the transversity distribution can be 
measured in Drell-Yan pair production with transversely polarized protons, where one 
combines the transversity quark and antiquark distributions in the two protons 0, 0, or 
in polarized jet production || |9|. In lepton-hadron scattering transversity is in principle 
accessible through semi-inclusive particle production, where the transversity distribution 
in the target is combined with a chiral-odd fragmentation function. In single-particle pro- 
duction one can measure the left-right asymmetries in the fragmentation of a transversely 
polarized quark, which are produced by T-odd fragmentation functions (Collins effect) 
fi~0| , [12|, 0, [14|, [16|]. In such processes there is usually a competing contribution in- 



volving a higher-twist distribution function in the target and a chiral-even fragmentation 
function, which makes extraction of the transversity distribution difficult. First estimates 
of the transversity distributions have been reported |T3L O based on the azimuthal asym- 
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metries measured by the HERMES |18[ and SMC [19j experiments. Alternatively, one 
can consider semi-inclusive production of two pions, where a T-odd structure in the frag- 
mentation function appears due to the possibility of interference between S- and P-waves 
in the hadronic final state |2(J. It has also been suggested to use so-called "handedness" 
correlations in multi-particle production to access the transversity distributions |Tj|. 

The transversity distributions have many interesting theoretical aspects, related to 
general properties of QCD as well as to the structure of the nucleon. In leading order 
(LO), certain inequalities have been derived by Soffer giving an upper bound on the 
transversity distribution in terms of the longitudinally polarized and unpolarized quark 
distributions; see also Refs. [ 22j . Their Q 2 -evolution and generalization to NLO have been 
discussed in Ref. 



3| Also, it has long been noted that in the non-relativistic quark model 
(no spin-orbit interaction) the transversity distributions are identical to the longitudinally 
polarized ones. In this sense, any measurements of the transversity distribution would 
have implications for our understanding of the relativistic structure of the nucleon. 

Recently a method has been formulated to calculate the quark- and antiquark distri- 



butions in the nucleon at a low normalization point in the large- N c limit |2J, |25| . In this 
limit the nucleon can be described as a soliton of an effective low-energy theory based 
on the dynamical breaking of chiral symmetry in QCD (chiral quark-soliton model) ||26|| . 
It has been shown that this fully field-theoretic description of the nucleon preserves all 
qualitative properties of the quark- and antiquark distributions known from QCD, such as 
positivity conditions, sum rules etc.. The unpolarized distributions have been computed 
in Refs. [24, 25, 2S]; the results are in good agreement with the well-known parameteri- 
zations of the parton distribution functions at a low input scale ||. In particular, this 
approach describes well the observed violation of the Gottfried sum rule and the flavor 
asymmetry of the unpolarized antiquark distribution |27|, [ft], [H], . The longitudi- 
nally polarized flavor-singlet distribution calculated in this approach |27) is also in good 
agreement with the parameterizations of the polarized DIS data ||. Interestingly, the 
large-iV c approach predicts a large flavor asymmetry of the polarized antiquark distribu- 
| |30], |3i"| , which could produce observable effects e.g. in semi-inclusive spin 
1021 



tion 24 



asymmetries 



and in Drell-Yan double-spin asymmetries with polarized protons 



the 



Also, the flavor-nonsinglet transversity distribution has been computed in Ref.f34 
isosinglet one in Ref. |27|, |2E| . We also note that the chiral quark-soliton model has been 
successfully applied to calculate skewed (non-diagonal) parton distributions [j35| . The 



model also describes well a large number of hadronic observables of the nucleon as well 
as the other octet and decuplet baryons, such as magnetic moments, electromagnetic and 
weak form factors, etc.; see Refs. |36|, |37| for a review. 

In this paper we report about a comprehensive investigation of the quark- and anti- 
quark transversity distributions in the large- N c approach of Refs. [^4], |25| . We calculate 
both the flavor-nonsinglet distribution, which appears in leading order of the 1/N C ex- 
pansion, and the flavor-singlet one, which is non-zero only in next-to- leading order. In 
Ref. |M| an approximation was used in dealing with the contribution of the Dirac contin- 
uum of quark states in the soliton; in this paper we perform a numerical calculation of 
the distribution functions without further approximations. This is particularly important 
for the antiquark transversity distribution. 

The effective chiral theory used to describe the nucleon is non-renormalizable and de- 
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fined with an explicit ultraviolet cutoff. An important problem when calculating parton 
distributions in this approach is to ensure that none of their essential properties such as 
positivity, sum rules etc. are violated by the ultraviolet regularization. For the unpolar- 
ized and polarized distributions this problem has been studied extensively in Refs. |24|, |25[] . 
Here we use the methods developed there to investigate the role of the UV cutoff on the 
transversity distributions. It turns out that these distributions are UV finite and can even- 
tually be computed without any cutoff. Still, in the actual calculation of the distributions 
one has to introduce a cutoff during the intermediate stages, and one must make sure that 
the finite result in the end does not depend on this intermediate regularization procedure. 
We explicitly address this problem here. In particular, we shall describe an interesting 
anomaly-type phenomenon observed in the calculation of the transversity distributions, 
which can be stated as the non-commutativity of the limits of infinite UV cutoff and the 
chiral limit (vanishing pion mass). A thorough understanding of this phenomenon is a 
prerequisite for a reliable calculation of the transversity distributions in the effective chiral 
theory. 

The plan of this paper is as follows. In Section |^ we summarize the definition and im- 
portant properties of the transversity distributions and the tensor charges within QCD. 
Section || gives a detailed exposition of the method used to calculate the transversity 
quark- and antiquark distributions at a low normalization point. After a brief descrip- 
tion of the effective low-energy theory and the chiral quark-soliton model of the nucleon 



in Section [D] we discuss in Section [T2] the 1/N C expansion of the transversity distribu- 
tions and the tensor charges. The calculation of the transversity distributions — both 
isovector and isoscalar - in the chiral quark-soliton model is described in Section \5.'d\ In 
particular, in this section we derive expressions for these distributions in the form of sums 
over contributions of quark single-particle levels in the soliton field which are used in the 
numerical calculations. Section [3^ deals with the important issue of ultraviolet regulariza- 
tion. The dependence of the transversity distributions on the ultraviolet cutoff is studied 
using analytic methods (gradient expansion). In Section pT5| we describe the anomaly-type 
phenomenon in the tensor charge, which is both of general theoretical interest as well as 
of importance in the numerical calculations. The nature of this phenomenon is clarified, 
and an explicit expression for the anomalous difference of the sums over contributions of 
occupied and non-occupied quark levels is derived. The corresponding difference for the 
full ^-dependent distribution is given in Appendix [A]. In Section [| we describe the results 
for the transversity distributions at the low scale. In Section we compare the results 
for the transversity distributions to those for the longitudinally polarized ones [ 3SJ, ^ 
at a qualitative level, considering various limiting cases of the chiral quark-soliton model 
which correspond to the non-relativistic quark model or the Skyrme model. Numerical 
results are presented in Section |4^ , where also a brief description of the numerical method 
used to evaluate the sums over quark levels is given. In Section O we demonstrate that 
the calculated transversity distributions satisfy the recently derived "large- N c versions" 
of positivity and Soffer inequalities ||40|| , which impose stronger constraints on the distri- 
butions than the "usual QCD" positivity and Soffer inequalities. Finally, in Section [5] we 
use our results for the transversity distributions to make predictions for the double spin 
asymmetries in Drell-Yan pair production in polarized pp and pp scattering. In particular, 
we compare the asymmetries calculated from our model distributions with those obtained 
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by assuming 5q(x) = Aq(x), an approximation frequently made in the literature. We 
find that our model distributions result in significant deviations from what is obtained 
with that approximation, which makes us hope that RHIC measurements may be able to 
discriminate between the two predictions. 

The transversity distributions at a low normalization point have been studied in a 
variety of other approaches. Model calculations of the distributions have been performed 
within the bag model M and a constituent quark model ET[. Furthermore, these distri- 



butions have been studied in the chiral quark-soliton model [27 , Q and a closely related 
description of the nucleon as a chiral soliton of the Nambu-Jona Lasinio model |^2| . The 
transversity distributions have also been estimated using QCD sum rule techniques |43|[ . 
For the tensor charges of the nucleon estimates from lattice simulations [|5[] and QCD 
sum rules [|B| are available. The tensor charges have also been calculated within the 
the chiral quark-soliton model |47j (the relation of these results to the present calcula- 
tion is discussed in Section £|7^). The evolution (scale dependence) of the transversity 
distributions has been studied by a number of groups, see Refs. HR |49|, |5(| . 



2 Transversity distributions in QCD 

In QCD, parton distributions are defined as expectation values of certain twist-2 light- 
ray operators in the nucleon state. The helicity decomposition of the matrix element of 
a general quark bilinear in a polarized nucleon ("spin density matrix") contains three 
twist-2 structures, which can be identified as the unpolarized, longitudinally polarized, 
and transversity distribution The explicit expression for the transversity distribution 
is 

<%(z) = l^e ixP+z -(N(P),S ± \M^h + lB^±[0,z]M^ \N(P),S ± ) Z+=Z±=0 , 

" (2-1) 
where ipf is the quark field of flavor f,z& light-like distance (z^z^ = 0), and the light-cone 
coordinates are defined as[| 

z ± = ^ = (2.2) 

Furthermore, [0, z] denotes the link operator (path-ordered exponential of the gauge field) 
required by gauge invariance. The nucleon state in Eq. (|2.1| ) is transversely polarized 
(Sj_ = —1, P^SI = 0), and £ ± = 7^Sj\ In the nucleon rest frame, the components of 
the polarization vector are given by 

Si = (0,0,0, A) (2.3) 



where A = 2Ss = ±1 and S3 is the spin projection on the 3-axis. In this frame Eq. (|2.1|) 
reduces to 

5g / (^) = -(253)/^e* (iV,53|^(0)(7 + 7 1 )757 3 [0^]^/(^)|iV,5 3 ) ^ . (2.4) 



=0 



x It is convenient here to define the light-cone vector components using the 1- rather than the 3- 
direction. With this choice the direction of transverse polarization of the nucleon can be chosen as the 



3-direction, which will make the expressions in Subsection 3.3 look more conventional 
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The operator in Eg. ( |2.1| ) is C-odd, so the corresponding antiquark distribution is 
obtained from the expression on the R.H.S. of Eq. fl2.1|) by 

8q f (x) = -5q f (-x). (2.5) 

Furthermore, the light-ray operator in Eq. ( |2.1|) is scale dependent, so that the distribu- 
tions depend on the normalization point. This dependence can be described by evolution 
equations for the transversity distributions, which have been studied in Refs. |48|, |49|, |50[ . 

The normalization integral for the distribution Eqs.( [2.1|) and ( [2.5|) is given by the 
so-called tensor charge of the nucleon, 

5qt = dx [5qf(x) — Sqt(x)] (2.6) 
J o 

which is defined as the matrix element of the local pseudo-tensor operator, 

(N{P),S\i>{-i)a^ l5 iP\N{P),S) = {P»S V -P v S»)5q f . (2.7) 
In particular, in the nucleon rest frame it is given by 

(N,SSl°l 3 l^\N,S 3 ) = 2M N (2S 3 )6q f . (2.8) 



3 Transversity distributions from the effective chiral 
theory 

3.1 The nucleon as a chiral soliton 

It is generally believed that in the large- N c limit QCD becomes equivalent to a theory of 
mesons, with baryons emerging as solitonic excitations While this connection alone 



has many interesting qualitative implications, it is not known at present how to derive this 
effective theory from QCD in full generality. However, quantitative calculations within 
the large-iV c ideology can be done in certain limiting cases where the dynamics of QCD 
simplifies. It is known that at energies far below the mass of mesonic resonances (say, the 
rho meson) the dynamics of strong interactions is governed by the spontaneous breaking 
of chiral symmetry. In fact, in the long wavelength limit the dynamics is completely 
described by the chiral Lagrangian containing only the Goldstone boson (pion) field to 
some given order in derivatives. To describe the nucleon as a chiral soliton, however, one 
needs an effective theory valid in a wider region of momenta. Such a theory has been 
derived within the framework of the instanton description of the QCD vacuum, which 
provides a "microscopic" picture of the dynamical breaking of chiral symmetry in QCD 
|52|j . In this approach the effective action for the pion field is obtained in the form of 



an integral over quark fields, which have obtained a dynamical mass in the spontaneous 
breaking of chiral symmetry, and which interact with the pion field in a chirally invariant 



way |52|, |53[: 



exp{iS cR [U(x)}) = / DipDip exp i / d^x^ - MU^ 5 )ip 



(3.1) 
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Here, are the fermion fields, M is the dynamical quark mass, and the Goldstone 
boson field is parameterized as a unitary matrix, U (x) , with 



W 5 {x) 



1+75 



U{x) 



1 - 75 



(3.2) 



In the long-wavelength limit, when expanding in derivatives of the pion field, the effective 
action Eq.( |3.1| ) reproduces the chiral Lagrangian with correct coefficients, including the 
Gasser-Leutwyler 0(p A ) terms and the Wess-Zumino term. However, the validity of the 
theory defined by Eq.( |3.1| ) is not restricted to the long-wavelength limit. 

In the effective theory derived from the instanton vacuum the dynamical quark mass 
in Eq. (|3.1|) is momentum-dependent and drops to zero for momenta of the order of the 
inverse average instanton size, ~ 600 MeV. This provides for a natural ultraviolet 
cutoff of the effective theory. In the calculation of quark distribution functions here we 
take the dynamical quark mass in Eq.( |3.1| ) to be constant and simulate the ultraviolet 
cutoff implied by the instanton vacuum by applying an external ultraviolet regularization 
to divergent quark loops. We shall show in detail that this is a legitimate approximation 
for the quantities considered here. 

The effective action allows to compute hadronic correlation functions at low energies 
in the large- N c limit. In particular, the nucleon in the effective theory at large-iV c limit 
characterized by a classical pion field which binds the quarks (chiral quark-soliton model) 
In the nucleon rest frame the classical pion field is of "hedgehog" form, 



26 



cxp 



-P(r) 



(3.3) 



where P(r) is called the profile function, with P(0) = — 7r and P(r) — > for r — > oo. 
The quarks are described by single-particle wave functions, which are the solutions of the 
Dirac equation in the background pion field, Eq . (|3~3|) : 



where 



H{U cl )\n) = E n \n), 



H{U cl ) = -2 7 V<9fc + M 7 °?7 



'75 

cl 



(3.4) 



(3.5) 



is the single-particle Dirac Hamiltonian in the classical background pion field. The spec- 
trum of H{U C \) includes a discrete bound-state level, whose energy is denoted by E\ ev , as 
well as the positive and negative Dirac continuum, polarized by the presence of the pion 
field. The soliton profile, P(r), is determined by minimizing the static energy of the pion 
field, which is given by the sum of the energy of the bound-state level and the aggregate 
energy of the negative Dirac continuum, the energy of the free Dirac continuum (U — 1) 
subtracted 
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Etot[U c 



E En- E e^ 



n 

occup. 



n 

occup. 



N c E lcv + N c ]T (E n -E®), 

n 

neg.cont. 



(3.6) 



and in the leading order of the l/iV c -expansion the nucleon mass is given simply by the 
value of the energy at the minimum, 



M N = E tot [U cl \. (3.7) 

The expression for the energy of the pion field, Eq. (|3.6|) , contains a logarithmic ultravi- 
olet divergence due to the contribution of the Dirac continuum and requires ultraviolet 
regularization (see below). 

In higher order of the l/iV c -expansion one must take into account the fluctuations of 
the pion field about its saddle-point value. A special role is played by the zero modes 
of the pion field. The minimum of the energy, Eq.( [3.6|) , is degenerate with respect to 
translations of the soliton field in space, and to rotations in ordinary and isospin space [for 
the hedgehog field, Eq. (|3~3|) , the two types of rotations are equivalent]. Quantizing these 
zero modes modes gives rise to nucleon states with definite momentum and spin/isospin 
quantum numbers ^4], |26|. This is done by subjecting the hedgehog field, Eq.( |3.3| ), to 



time-dependent translations and flavor rotations, 

f/ cl (x) -> R(t)U cl (x-X(t))rf(t), (3.8) 

where R(t) is an SU(2) rotation matrix, and computing the functional integral over the 
collective coordinates within the l/iV c -expansion. The functional integral over transla- 
tions, X(i), can be reduced to a Hamiltonian system describing the free motion of the 
soliton center-of-mass, with mass Mjy, Eq.( |3.7| ), whose eigenfunctions are plane waves 
with given three-momentum, exp(zP-X). [The 0(1/N C ) contribution of the collective 
translations to the nucleon mass can be neglected in the applications considered here.] In 
the functional integral over collective rotations, R(t), the rotational energy is given by 

-^rotating soliton -^static soliton — I tl [f2(t)f2(t)] + . . . , (^-9) 

where 

n(t) = fi a (*)y = -iR\t)R(t) (3.10) 

is the angular velocity, and / is the moment of inertia of the soliton, which is given by a 
double sum over quark single-particle levels in the background pion field, 

T N c ^ ^ (n\r a \m)(m\r a \n) 

" n m ■l- J m - C/ n 

occup. non— occup. 

Here the sum over n runs over all occupied states, i.e., the discrete level and the negative 
Dirac continuum, the sum over m over all non-occupied states, i.e., the positive Dirac 
continuum. (The ultraviolet regularization of this quantity will be discussed below.) It is 
important that the moment of inertia is of order N c , so the typical angular velocities are 

m - y (3.i2) 



9 



and one can compute the functional integral by expanding in powers of the angular ve- 
locity. To leading order in 1/N C the collective motion is described by a Hamiltonian 



C2 rp2 

H„, = £ = f v (3.13) 

where S a and T a are the right and left angular momenta, and the Hamiltonian Eq. ( |3.13| ) 
has been obtained by the "quantization rule" 

fia -> y- (3.14) 

This Hamiltonian describes a spherical top in spin/isospin space, subject to the constraint 
S 2 = T 2 , which is a consequence of the "hedgehog" symmetry of the static pion field, 
Eq.(|3.3|). Its eigenfunctions, classified by S 2 = T 2 ,S 3 and T 3 are given by the Wigner 



finite-rotation matrices 26 



4%£(R) = V2STT(-ir^D s _^ S3 (R). (3.15) 

The four nucleon states have S = T = 1/2, with 5*3, T 3 = ±1/2, while for S = T = 3/2 
one obtains the 16 states of the A resonance. The rotational energy, S(S + 1)/(2J), gives 
a l/iV c -correction to the nucleon mass, which should be added to Eq. (|3.7| ). In particular, 
the nucleon-A mass splitting is given by 

M A -M N = A (3.16) 



The saddle-point solution, Eq.(|3.3j), and the collective quantization procedure outlined 
above, not only give rise to nucleon states of definite momentum and spin/isospin quantum 
numbers in the effective chiral theory defined by Eq. fl3.lp , they also imply a prescription 
for the calculation of matrix elements of arbitrary composite quark operators between 
nucleon states. For a review of the applications of this model to baryon observables such 
as masses, form factors, etc. we refer to Ref . p6| . We shall apply this approach to calculate 
the nucleon's transversity distributions in Subsection |3.3| . 

3.2 The transversity distributions in the large— N c limit 

Before embarking on the calculation of the transversity distributions in the chiral quark 
soliton model it is useful to establish the large-iV c behavior of these distributions on 



general grounds 0. Standard N c counting tells us that at large iV c the tensor 



charges of the nucleon scale as 

8u-8d ~ iV c , 8u + 8d ~ 1, (3.17) 

i.e., the isovector matrix element is leading relative to the isosinglet one. This behavior 
is analogous to that of the axial charges, which scale as 

gf = Au-Ad ~ N c , gf ] = Au + Ad ~ 1. (3.18) 
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Combining Eq. fl3.17|) with the fact that in the large-iV c limit the parton distributions 
are concentrated at values of x of the order of 1/N C one obtains that the transversity 
distributions at large N c scale as 

5u(x)-8d(x), 8u(x)-5d(x) ~ N 2 C f(N c x), (3.19) 
5u(x)+8d(x), 5u(x)+5d(x) ~ N c f(N c x), (3.20) 

where f(y) is a stable function in the large iV c -liniit, which depends on the particular 
distribution considered. 



3.3 Calculation of the transversity distributions 



We now outline the calculation of the transversity quark and antiquark distributions in 
the chiral quark soliton model. The methods for computing parton distributions at a low 
normalization point in this approach have been developed in Refs. p4|, p5| , |29"|| , and we 
refer to these papers for all general methodological questions. 

When computing quark/antiquark distribution in the effective chiral theory it is as- 
sumed that the normalization point of the distributions is of the order of the ultraviolet 
cutoff of the effective theory, i.e. of (9(600 MeV). At this scale the QCD twist-2 quark 
operators can be identified with the corresponding operators in the effective theory. The 
gluon distributions are zero at this level of approximation [2~3|. A more explicit justifica- 
tion for this procedure is provided by the instanton picture of the QCD vacuum, which 
allows to derive the low-energy effective theory. It was shown that in leading order of the 
packing fraction of the instanton medium the quark and antiquark distributions saturate 
the nucleon momentum and spin sum rule, while the gluon distributions are zero |2"3| , . 
The calculations of parton distributions in Refs. [0, |25|, |56], [38|, |39| were performed at 
this level of approximation. In the case of the transversity distributions, since there is no 
transversity gluon distribution, one may expect this "quarks-antiquarks only" approxi- 
mation to give even better results than in the unpolarized and longitudinally polarized 
case. 



Having expressed the QCD operators in Eqs.( |2.4j , |2.5| ) in terms of the quark fields of 
the effective theory, we can now compute their nucleon matrix elements within the 1/N C - 
expansion, using standard techniques. It is convenient to work in the nucleon rest frame, 
where the classical pion field describing the nucleon (up to collective translations and 
rotations) is given by Eq.( |3.3| ); in this frame the matrix elements defining the transversity 
distributions take the form Eqs.( [2.4| , |2.5|) . Matrix elements of quark bilinear operators 
such as Eqs 



2Tq ) can be reduced to those of time-ordered products of quark fields, 
which can be calculated with the help of the Feynman Green function of the quarks in 
the background pion field, 

GV(2/ ,y;x°,x) = (^y|[^-tf(l0rV,x>. (3.21) 

Here the saddle-point pion field is the slowly rotating hedgehog field, Eq. ( p.8[ ). For this 
ansatz the Green function, Eq. (|3.21[ 



takes the form 



[id t -H{U)\ 



R(t) {idt-HiU^-ttit)]- 1 B}{t), 



(3.22) 
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where Q(t) is the angular velocity, Eq.( |3.10D . Performing the functional integral over 



collective coordinates of the saddle-point field as described in Subsection pTT] , projecting 
on nucleon states with definite momentum and spin/isospin quantum numbers one obtains 
the following "master formula" for the expectation value of a color-singlet time-ordered 
quark bilinear operator in the nucleon: 

(P = 0, S = T, S 3 , T 3 \ T {^(x)T^(y)} |P = 0, S = T, S 3 , T 3 ) 

R(T)=R 2 

k T=S 



2M N /cf X dRi dR 2 



(R2 



x(-i)N c Tr 



R\x°)TR(y ) (y°,y-X 



^K^i) J DRDet[idt-H(U d )-n(t)] 

R(-T)=Ri 
1 



%d t - H(u cl ) - n(t) 



|z°,x-X) 



(3.23) 



Here T denotes a matrix in Dirac spinor and isospin space, and Tr . . . implies the trace 
over Dirac and flavor indices (the sum over color indices has already been performed). 

The functional integral over rotations, R(t), in Eq. fl3.23|) can be computed by ex- 
panding the integrand in the angular velocity, Q, Eq.( p.l0|) , which is of order 1/N C . The 
expansion of the determinant of the Dirac operator gives rise to the "kinetic term" 



Det [id t - H(U c i) - Q(t)] oc exp 



Y jdt nl(t) + 



(3.24) 



where I is the moment of inertia of the soliton, Eq.( [3.1l|) . The functional integral with 
this action can now be computed exactly; it corresponds to a rigid rotator described by 
the Hamiltonian Eq. (|3.13|) . In addition, one has to expand the quark Green function in 
Eq. (|3.23|) in powers of Q. The minimum power of Q required to obtain a non-zero result 
determines the order of the matrix element in the 1/N C expansion, which is in general 
different for the different isospin components of the matrix element of a given operator. 

Isovector transversity distribution. We now apply the above prescription to the calcu- 
lation of the transversity distributions, Eqs.( [2.4| , [2.5|) . It turns out that in the case of the 
isovector transversity distribution the R.H.S. of Eq. ( |3.23| ) is non-zero already in zeroth 
order of the expansion of the quark propagator in Q(t), in agreement with the fact that 
this isospin component is leading in the l/iV c -expansion. The functional integral over 
rotations in leading order of 1/N C (i.e., neglecting all time dependence of the collective 
rotation) just produces a delta function which enforces Ri = R 2 . The result can be 
written in the form 

8u(x) - 8d(x) = -(2S 3 )2M N J d 3 X J dR (j)* T3 s 3 (R) /=1 (X, R;x) ^ Sa (i2)(3.25) 

where I=1 is an operator in the space of functions of the collective coordinates 

„o_,i 



O i=1 (X,R;x) 



R^t 3 R i±^ 757 3 



X (2°,Z-X| 



idt - H{U ( 



clJ 



|0,-X> 



(3.26) 



-2°, Zj_=0 
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The integral over rotations is readily performed; one introduces the rotation matrix in 
the adjoint representation, 



tfr 3 R = D 3a (R)r a , 
and makes use of the identity 



D ba {R) = ^Ti{r b Rr a R^, (3.27) 



dR^*(R)D 3a (R)^ S3 T 3 (R) 



--(2T 3 )(2S 3 )5 a3 . 



(3.28) 



For further evaluation one writes the quark Green function in the static classical pion 
field in Eq.( |3.26|) in frequency representation, 



XI 



id t - H(U cl ) 



I"- -X. ) / ^ e— °(z - X| | - X ) (3.29) 

' 2lX oo — H {U c \) 



(the choice of contour for the u;-integral will be discussed below), and brings the matrix 
element into diagonal form by introducing the finite-translation operator, 



(z-X|...| - X) 



(-X| exp(iP k z 



X), 



(3.30) 



where P k denotes the momentum operator in the space of quark single-particle wave 
functions, Eq. fl3.4|). In this way one obtains 



5u(x) — 5d(x) 
d 3 XTr 



(2T 3 ) {-%) 



duj 
2^ 



x 



r 3 l+ J° 7l 7 5 7 3 (-X| 5{xM N - uj - P 1 ) 



H(U ( 



X) 



clJ 



.(3.31) 



The delta function here is the result of integrating the exponential factors in Eqs. (|3.26 ), 
Q3.29| ) and ( ^.30[ ) over z°, keeping in mind the constraint z 1 = —z°. In fact, the R.H.S. of 
Eq. (|3.31|) has the form of a functional trace of an operator in the space of quark single- 
particle states, and can be written more concisely as (Sp denotes the functional trace) 



Su(x) — 8d(x) 



(2T 3 ) (-«) 

t + tV 



did 

2^ 



xSp 



T 



l5l 3 5(xM N -id-P 1 )- 



id 



H(U C] 



(3.32) 



In the above expressions vacuum subtraction is implied, i.e., one should subtract the 
corresponding expressions in which the Hamiltonian in the hedgehog pion field, Eq.( |3.3|) , 
is replaced by the free Hamiltonian (U = 1). 

Eqs.( |3.3l| ) and ( |3.32j ) serve as the basis for the actual computation of the isovector 
distribution. An explicit expression, suitable for numerical calculations, can be derived 
by substituting in Eq.( |3.32j ) the spectral representation of the quark Green function in 
the background pion field, 



uo - H(U C 



E 



occup. 



|n)(n| 



E n -i0 



+ 



id 



E 



|n)(n| 



id — E n + iO 



(3.33) 



non— occup. 
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Here, E n and \n) are the single-particle eigenvalues and eigenstates of Eq. ( |3.4j ). The 
poles are shifted to the upper /lower half of the u;-plane, corresponding to whether the 
single-particle levels are occupied or not. Substituting Eq.( ^.33[ ) into Eq.( |3.32j ) we obtain 
a representation of the isovector distribution as a sum over quark single particle levels in 
the background pion field: 



8u{x)-8d{x) = (2T 3 )^^ 



/ 



E - E 



V.- 



v occup. non— occup. , 

x(n|r 3 l+ ^° 7l 757 3 5(xM N - E n - P l )\n). (3.34) 

For many purposes it is desirable to have representations of the distribution function 
involving sums over only occupied, or only non-occupied, levels. Such representations can 
be obtained from Eq.( |3.34! ) if one notes that the sum over all levels of the matrix element 
on the R.H.S. of Eq.(|3~34T) is zero: 



EH^ 1 + 2° 7l ^7 3 S(xM N - E n - P 1 )^) = 0. (3.35) 

n all ~~ 

This condition is actually equivalent to the locality condition of the anticommutator of 
quark fields in the effective low-energy theory, which ensures that one gets the same 
result for the distribution function in the chiral quark-soliton model if one starts from 
the QCD definition as the matrix element of the quark bilinear ip(0) . . .ip(z), Eq.( [2.1|) , or 
from the equivalent definition where as the matrix element of the QCD operator where ip 



and ip have been anticommuted, — ip(z) . . . |24], |25 |. It is crucial that the ultraviolet 
regularization of the effective low-energy theory does not destroy the property Eq.( |3.35| ), 
as has been discussed extensively in the context of the isoscalar unpolarized and isovector 



polarized distributions in Refs.El, A detailed investigation of the conditions under 



which Eq. (|3.35|) holds in the case of the isovector transversity distribution is presented in 
Subsection |3.5| and Appendix here we shall simply take this property for granted. In 
particular, Eq. (|3.35|) can be read as saying that 

(n\ . . . \n) = — (n\...\n). (3.36) 

n n 
occup. non— occup. 



This allows to write instead of Eq.( [3.34 ) 



Su(x) — 5d(x) 



(2T3) 2A^v £ {n ^l±p^^ f f S(xMN _ En _pi )ln) (3. 37) 

On. ^ 



occup. 

,,0^,1 



-(2T 3 )^^ E (n\r 3l -±^^5(xM N -E n -P')\n). (3.38) 



3 

non— occup. 



All three representations of the isovector transversity distribution, Eqs. (|3.34|) , (|3.37|) and 



( |3.38| ) are equivalent provided Eq. ( f3.35| ) holds. In all cases, the corresponding expressions 
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for the antiquark distributions are obtained from those for the quark distributions by the 
substitution Eq. (|2.5| ). It is understood that one should subtract from Eqs. fl3.34|) , (|3 . 3 7| ) 
and ( |3.3<j| ) the corresponding sums over eigenstates of the free Hamiltonian (U — 1). 

The expressions Eq. (|3.37|) and ( |3.38| ) can be used for numerical evaluation of the isovec- 
tor distribution function, see Section [4.2| . In particular, we shall verify the equivalence of 
the two representations in the numerical calculations. 

Isosinglet transversity distribution. The calculation of the isosinglet transversity dis- 
tribution is slightly more complicated than that of the isovector one. In the isosinglet case 
a non-zero result is obtained only after expanding the integrand in Eq. fl3.23|) to first order 
in the angular velocity, Q. As a consequence this distribution is suppressed by a factor 
1/N C relative to the isosinglet one, in agreement with the general iV c -counting arguments 
of Subsection \J72. 



Quark distribution functions are given by matrix elements of non-local operators 
(quark bilinears at different space-time points). The ^-expansion of matrix elements 
of such operators poses some special problems, which have been discussed extensively 
in Ref.[p9[1 in connection with the calculation of the isovector unpolarized distribution. 
Terms of first order in Q(t) arise from the first-order expansion of the Green function in 
Eq. (|3.23|) , as well as from the expansion of the structure 



R\x°)TR(y ), 

which is non-local in time.0 There are thus two types of contributions to the distribution 
functions. For the isosinglet transversity distribution the result obtained after expanding 
to first order in Q can be written as [cf. Eq. (|3.25| )1 



Su(x) + 5d(x) 



(2S 3 )2M N J d 3 X J dR <P* T:iS3 (R) 

+ /=o, (2)1 (x,R,S;x)<f> T3S3 (R), 



(3.39) 



where 



O I=0 'W{X,R,S;x) 



x (z°,z-X| 



— S a T a 



id t -H(U cl ) V2J ■ J id t -H(U ch 

dz° ... 



|0,-X) 



(3.40) 



J z 1 =-z°, z±=0 



O I=0 ' ( 2 )(X, R, S; x) = (-i)N c [ ^e ixM » z ° z° Tr 
1 



(3.41) 



x (z u ,z-X| 



id t - H(U ( 



|0,-X) 



ci; 



-z°, z ± =0 



Now the operators acting on the wave functions in collective coordinates involve also 
the spin operator, S, which arises from replacing the angular velocity according to the 
"quantization rule" , Eq. ( |3.14| ) . 

2 If one computed not the distribution function directly but its moments, the second type of contribu- 
tion would arise from the presence of derivatives acting on the quark fields in the local twist-2 operators, 
which in the chiral quark-soliton model become time derivatives acting on the rotational matrices in 
Eq.(p3|). 
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Note that in Eqs.( |3.40D and ( |3.41| ) the spin operators do not commute with the ro- 
tational matrices, so in principle one should be careful about their ordering. However, 
it turns out that in the case of interest here there is no ordering ambiguity. By explicit 



calculation, using the methods described in Ref.|29|], one can show that the commuta- 
tor terms corresponding to the differences between different operator orderings give zero 
contribution in the final result .0 

The integrals over rotational matrices in Eqs.( 3.40| ) and (3.41) are nothing but the 
average of the spin operator in the rotational state, 



J dR<j>* SsTa {R) S a (Ps 3 t 3 (R) = S 3 6 a3 . 



(3.42) 



The further evaluation of the expressions proceeds largely in analogy to the isovector case. 
Passing to the frequency representation of the quark Green functions, cf. Eq. (|3.29|) , and 
making use of Eq. Q3.3UD , the contributions Eqs.( |3.40| ) and ( |3.41| ) become 



[5u(x) + 5d(x)] {1) 
'1 + 7V 



N r M 



N 



41 



2^ 



xSp 



W 3 S(xM N -uj-P 1 ) 



co - H(U cl ) 



UJ 



H(U C 



[Su(x) + 5d(x)] {2) 



41 [ l) dxJ 2tt 
l5l 3 5(xM N - uj - P 1 )- 



UJ 



H(U C 



(3.43) 



(3.44) 



The derivative in x in Eq.( |3.44| ) results from the factor of z° present in Eq.( |3.41| ). One 
notes that this contribution is, up to a factor, equal to the derivative in x of the leading- 
order O(Q ) result for the isovector transversity distribution, Eq.( |3.32P , 



[8u(x) +5d{x)] {2) 







2IM N dx 



[Su(x) — 5d(x)] 



leading 



(3.45) 



Substituting here Eqs.( |3.37P or ( |3.38|) one obtains a representation of this contribution 
as a simple sum over (occupied or non-occupied) quark single-particle levels. For the 
contribution (1), Eq. (|3.43|) , a representation as a double sum over quark levels can be 



derived, following the steps outlined in Appendix A of Ref . p9f . Again one substitutes the 
spectral representations, Eq.( |3.33| ), for the two quark propagators. A subtle point is the 
occurrence of a double pole in the w-integral for those terms where the energies in the 
two denominators coincide, see Ref.f2T| for details. Assuming a quasi-discrete spectrum 
of levels, as appropriate for numerical calculations using a finite box, and separating 



3 The absence of ordering ambiguities can be shown to be a general feature when computing 1/N C - 
suppressed quantities, such as the isosinglet transversity distribution considered here, in leading nonva- 
nishing order, that is, at level O . Ordering ambiguities may arise, however, when considering 1/N C - 
corrections to quantities which are non-zero already in leading order of the l/A^-expansion, e.g. in 
f^-contributions to the isovector transversity distribution or to the isovector longitudinally polarized 
distribution. These difficulties cause a violation of the PCAC relation and are the subject of on-going 
investigations. 
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explicitly the contributions of levels with E m ^ E n and E m = E n in the double sum over 
levels n and m one obtains a representation of the isosinglet transversity distribution in 
the form 



[5u(x) + 5d(x)] {1) 



2 E 



N c Mn 
41 



( 



1 



E - E 

n n 
occup. non— occup. 

cu,i 



E — E 



(n\r 3 \m)(m\ 1 + 1 7 757 3 8{E n + P 1 - xM N )\n) 



V] (n\r 3 \m)(m 

m 
E —F, 



0„,1 



1 + 7 U 7 



l57 3 S'(E n + P :i -xM N )\n) 



> . 



(3.46) 



As in the case of the isovector distribution, one can show that the sum over all levels n 
of the terms in braces in Eq.( |3.46| ) is zero (see Subsection |375| and Appendix [A]), 



£{•••} 



o. 



(3.47) 



all 



so that we can write Eq. (|3.46|) equivalently as 



[5u{x) + 5d{x)} {1) 



N c Mn 
21 

N C M N 



21 



E {•••} 

n 
ccup. 

E {•••} 



(3.48) 
(3.49) 



non— occup. 



The expressions Eqs.( |3.48"D and ( |3.49| ), as well as Eq.( |3.45|) together with Eqs.( |3.37|) and 
( p.38|) , will be used in the numerical calculation of the isosinglet transversity distribution 



in Section 472 



3.4 Ultraviolet regular izat ion 

The effective chiral theory used to describe the nucleon is non-renormalizable and under- 
stood to be defined with an explicit ultraviolet cutoff. When calculating parton distri- 
butions in this approach one must ensure that none of their essential properties, such as 
positivity, sum rules etc., are violated by the ultraviolet regularization. For the isoscalar 
unpolarized and isovector polarized distributions this problem has been studied in detail 
in Refs. f24], p5f| . These distributions were found to contain ultraviolet divergences which 
require regularization, and it was shown that a possible regularization, preserving all 
important properties of the parton distributions, is by way of a Pauli-Villars subtraction. 

We now discuss the issue of ultraviolet regularization in the case of the transversity 
distributions. Our investigation consists of two parts. First, we show that the expressions 
for the transversity distributions in the effective low-energy theory are in fact ultraviolet 
finite, and thus can in principle be computed without an ultraviolet cutoff. Second, we 
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investigate the role of ultraviolet regularization in the "locality" conditions, Eqs.( |3.35| ) and 
( |3.47 )) which ensure the equivalence of the representations of the transversity distributions 
as sums over occupied and non-occupied quark levels. In Subsection [J75] we describe an 
interesting anomaly-type phenomenon found in these sums over quark levels, which takes 
the form of non-commutativity of the chiral limit and the limit of infinite ultraviolet 
cutoff. 

In order to investigate the dependence of the transversity distributions calculated in 
the effective low-energy theory on the ultraviolet cutoff, we consider the formal limit of 
large soliton size, in which one can derive explicit expressions for the distributions as 
functional of the classical pion field, U c \(x) (gradient expansion) |24], [25|]. The gradi- 
ent expansion can immediately be derived from the representations of the distribution 
functions as functional traces with the quark Green function, Eq. fl3.32|) in the isovector, 
and Eqs. (|3.43|) and (|3.44j ) in the isoscalar case, if one substitutes an approximate form 
of the quark Green function, appropriate for small gradients of the classical pion field, 
V fc tfd(x) <M (fc = 1,2,3): 



uo + H{U< 



cl 



1 



u - H(U cl ) tu' 2 
1 



cl) 



u 2 + V 2 -M 2 -iM(fU cl 



TV 



( W - l7 VV fc +M 7 ^) 



uj 2 +V 2 -M 2 



E 

n=0 



iM(fUZ 



lu 2 +V 2 -M 2 



(oj-i^ k V k + M 7 °[^ 5 ) . (3.50) 



In this formal expansion increasing numbers of gradients of C/ C i(x) come with increasing 
inverse powers of u and momentum, so that the leading ultraviolet divergences (if any) 
of the distributions can be read off from the leading-order gradient expansion. 

Isovector transversity distribution. In the case of the isovector transversity distribution 
the first non- vanishing contribution in Eq. (|3.32 ) comes from the term with n = 2 in the 
expansion Eq. (|3.50|) . Computing the resulting functional trace using a basis of plane-wave 
states one obtains the leading-order gradient expansion in the form 



[Su(x) — 5d(x)} 



grad. exp. 



d£ e 



/ /=1 (0, 



/ J=1 (£) 



N C M N M 



1 a 



X 



da jd(3 /^tr fl {f/ cl (z-<e 3 ) [^(z-/^)] [djU^z)} r k } 

(3.51) 



48 7T 3 



(e i2j 5 kl +e ll] 5 k2 ) + i(e lij 5 kl -e i2j 5 k2 ) 



We note that in deriving this result we have assumed that the classical pion field drops 
faster than 1/r 2 for r — > oo. This behavior is required in order to be able to drop certain 
surface terms involving the pion field at r = oo, which arise from terms in Eq.( |3.50| ) 
with n = and 1. Such contributions are at the heart of the anomaly-type phenomenon 
described in Subsection 3~ 5 and Appendix |A|, and we shall return to this point there. 

The gradient expansion, Eq.( |3.51| ), tells us that the isovector transversity distribution 
is ultraviolet-finite: the cu-integral in Eq. (|3.32|) is convergent. Thus this distribution does 
not require the ultraviolet cutoff of the effective theory and can be computed in the limit 
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of infinite cutoff. We further note that the gradient expansion result for the isovector 
distribution is real. By explicit calculation one can show that of f I=1 (£), Eq.( p.5l| ) the 
real part is even, the imaginary part odd in £. 

Isosinglet transversity distribution. In a similar way one can derive the gradient ex- 
pansion for the isosinglet transversity distribution. (We do not quote the expressions 
here.) The gradient expansion shows that this distribution is finite in the limit of large 
UV cutoff and does not require regularization. 



3.5 Anomaly— type phenomenon in the tensor charge 

The tensor charges, and, more generally, the transversity distributions show an interest- 
ing anomaly-type phenomenon, which we shall discuss now. Aside of being of general 
interest, this anomaly has direct implications for our procedure of calculating the transver- 
sity distributions in the chiral quark-soliton model, as it is related to the equivalence of 
the representations of the distributions as sums over occupied or non-occupied quark 
levels, Eqs.( [3.37D and (|3.38| ). This is equivalent to the condition that the sum of the 



single-particle matrix elements over all levels (i.e., occupied and non-occupied) be zero, 
Eq. (|3.35|) . A thorough understanding of this phenomenon is thus necessary for a reliable 
calculation of the distributions. In this section we illustrate the essential points by con- 
sidering the simplest case of the isovector tensor charge; the corresponding calculation for 
the x-dependent distribution function presents only technical difficulties and is presented 
in Appendix |A]. 

Anomaly in the tensor charge. At its most general, the anomaly we are dealing with 
is a statement about the functional trace of a chirally odd operator in the space of quark 
single-particle wave functions in a chiral background field. Let us consider the following 
functional trace: 

Sp [r a 7o i(J 0k 7s] , (3.52) 

which appears in the expression for the tensor charge. Naively this trace would be zero 
because the corresponding matrix traces in flavor and spin spaces are zero. However, 
actually we deal here with an uncertainty of the type • oo. To resolve this uncertainty 
one has to introduce a regularization. We choose here for illustrative purposes to regularize 
the trace ( |3.52|) in the following way: 

Sp [t°7o ia ok 7 5 ] reg = JimSp r a 7 ia ok 1^ H ' 2 . (3.53) 

Here 

H 2 (U) = -V 2 + M 2 +iM(fU 75 ), (3.54) 

is the Hamiltonian squared of the chiral quark-soliton model. In the limit of small e 
the expression in the RHS of eq. (|3.53|) can be computed analytically with the help of 
semiclassical expansion. The leading contribution has the following form: 

limSp [r a 7o ia ok j^ 2 } = / d 3 x d k tr [r a U(x)] . (3.55) 



We see that the trace ( |3.53| ) is linearly divergent provided the chiral field U{x) falls off 
not too rapidly at spatial infinity. For the soliton solution in the chiral limit the chiral 
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field behaves at the spatial infinity as U(x) — 1 ~ ^£ what leads to a non-zero coefficient 
in front of the linear divergence 1 / ^ft. The corresponding coefficient can be computed in 
terms of axial charge of the nucleon keeping in mind the asymptotics of the soliton 
solution in the chiral limit: 



8ttF2 x 



The result has the form: 



lim Sp r a 7 o «r * = K . (3.57) 



Note however that in the case of non-zero quark masses (whatever small they are) the 
linear divergence and hence the anomaly is zero. For the calculations in the finite volume 
the condition for the absence of the anomaly is M^xbox size^> 1. When computing the 
tensor charge (or the transversity distribution) by summing over quark levels in a finite 
volume, in the strict chiral limit one would find that the summation over occupied levels, 
Eq. (|3.37|) , and non-occupied levels, Eq.( |3.38|) , does not give equivalent results; rather, the 



difference diverges linearly with the ultraviolet cutoff. This is indeed what we observe in 
the numerical calculations. 

Elimination of the anomaly in the numerical calculations. In the numerical calcula- 
tions, which are based on diagonalization of the Dirac Hamiltonian in a finite box, it 
is important to eliminate the "anomaly", Eq. fl3.57p . A way to do this is to modify the 
soliton profile at large distances in a way that it vanishes faster than 1/r 2 . For example, 
one may use a modified version of the variational profile suggested in Ref. | 



PmA t ) = _2 arctan 



2 

(3.58) 



pf(l + M w r) exp(-M w r) 



where tq « l.OM" 1 is the usual soliton size parameter. For M n ^ this profile decays 
exponentially at large r, while for small r it differs only very little from the massless 
profile. The "anomaly", Eq. (|3.57| ), can then be made arbitrarily small by choosing the 
radius of the box used in the numerical calculation, D, sufficiently large so that 

M W D > I. (3.59) 

The limit M n — ^ is then taken at the very end of the calculation, by extrapolation 
of the numerical data. A measure of the successful elimination of the "anomaly" in the 
numerical calculations is the equivalence of the results obtained summing over occupied 
and non-occupied states. In the actual calculations we have used values of M T of the 
order of 0.3. . . 0.6 M, and a box radius of 20M -1 . For these parameters we observed 
equivalence of summing over occupied and non-occupied states at the level of 1 . . . 2%. 



4 Results 

4.1 Transversity vs. longitudinally polarized distributions 

Having derived the expressions for the transversity distributions in the effective low- 
energy theory, we now proceed to compute the distributions and discuss their properties. 
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Before turning to the numerical evaluation of the expressions it is instructive to compare 
the expressions for the transversity quark- and antiquark distributions with those of the 
corresponding longitudinally polarized distributions at a qualitative level. In Section W72 
we shall then compare the numerical results. 

Since the nucleon's axial and tensor charges are the same in the non-relativistic quark 
model, it is generally thought that differences between the two distributions are a measure 
of "relativistic effects" in the nucleon. The chiral-quark soliton model, which is a fully 
relativistic description of the nucleon, offers the unique possibility to study "relativistic 
effects" in the quark- and antiquark distributions in a controlled way. The relevant 
parameter here is the "soliton size", i.e., the radius of the classical pion field of the 
nucleon. Although in reality the soliton size is determined by the minimization of the 
classical energy, Eq.( |3.6| ), it is instructive to regard it as a parameter and to study the 
dependence of quantities such as the quark and antiquark distribution functions on it |58 



For small soliton sizes the bound-state level of quarks becomes weakly bound, and the 
lower Dirac components in its wave function become small. It was shown that in this limit 
nucleon matrix elements of a variety of local operators tend to their corresponding values 
in the non-relativistic quark model [[58 1. On the other hand, in the limit of large soliton 
sizes the bound state level approaches the negative continuum, and one may perform an 
expansion of nucleon matrix elements in inverse soliton size, which technically is obtained 
by expanding in gradients of the classical pion field (gradient expansion). In this limit 
this picture of the nucleon shows many similarities with the Skyrme model. In this sense 
one may say that the chiral quark-soliton model interpolates between the non-relativistic 
quark model and a Skyrme soliton picture of the nucleon. 

The quark model limit (small soliton size). It is thus interesting to study the transver- 
sity and longitudinally polarized distributions in dependence on the soliton size. Consider 
first the isovector distributions, which are leading in the l/iV c -expansion and both given 
by simple sums over quark single-particle levels, see Eq.( [3.37|) and Ref.|[24||. For small 
soliton size the contributions due to the polarized Dirac continuum become negligible, 
and both sums are dominated by the contribution from of the bound-state level: 
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Evaluating this contribution using the explicit form of the bound-state level wave function 
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(4.3) 



where and j(k) are the Fourier transforms of the radial wave functions corresponding 
to the upper and lower components of the Dirac spinor wave function (see Appendix B 



of Ref. |p4| for details). From the Dirac equation in the hedgehog pion field it follows that 
in the limit of small soliton size the lower Dirac component of the level wave function 
becomes small: j(k) — > 0, so that the bound state effectively becomes non-relativistic. 
One sees that in this case the expressions for the isovector transversity and longitudinally 
polarized distributions coincide. 

In a similar way one may investigate the limit of small soliton size in the isoscalar 
polarized distributions. This is slightly more complicated, as these distributions are given 
by double sums over quark levels. Consider the representation of the isoscalar distribution 
in the form of a spectral integral, Eqs.( |3.43"D and Eq. (|3.44|) . One can show that in the limit 
of small soliton size the dominant contribution is the sum of two contributions which one 
obtains by replacing one of the propagators in Eqs. (|3.43| ) by the pole associated with the 
bound-state level, |lev)(lev|/(o; — E\ ev ), the other one by the free propagator, l/(o> — M). 
In this way one can show that also in the isoscalar case the transversity and longitudinally 
polarized contributions tend to the same function in the "quark model" limit. 

The "shyrmion" limit (large soliton size). In the formal limit of large soliton size the 
bound-state level of quarks disappears in the negative-energy Dirac continuum. In this 
limit functions of the soliton field can be computed by expanding in gradients of the pion 
field; more precisely in the parameter dU/M. The techniques for deriving the gradient 
expansion of quark and antiquark distributions have been developed in Refs. |25j . For 
the isovector transversity distribution the result of the leading-order gradient expansion 
has been given in Eq.( |3.51 ). The corresponding expression for the isovector longitudinally 
polarized distribution has been derived in Ref. |24 
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(In the last line we have given the result in the case when the logarithmic divergence of 
the integral is regularized by a Pauli-Villars subtraction.) Comparing the two expressions 
we see that in the limit of ro — > oo (tq is the parameter characterizing the soliton size) the 
isovector transversity distribution, Eq.( |3.51| ), is suppressed relative to the longitudinally 
polarized one, Eq.( [4.4D , by a factor of l/(Mr ). Thus we see that in this limit the two 
distributions behave differently at a qualitative level. This is not unexpected, since, as 



4 The isovector longitudinally polarized quark- and antiquark distributions coincide only in the hypo- 
thetical limit of large soliton size. For finite soliton sizes they are, of course, different. 
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explained above, the limit of large soliton size can be regarded as the opposite of the 
non-relativistic limit. 

In the previous calculation of the isovector polarized quark distributions in the chiral 
quark-soliton model p4[ it was observed that the gradient expansion expression, 
when evaluated at the the physical soliton size, gives a fairly realistic description of the 
isovector polarized antiquark distribution. (The quark distribution, in contrast, is poorly 
described by the gradient expansion expression.) The same can be expected to apply 
to the isovector transversity distribution. If this is so, it implies that we should expect 
the isovector transversity antiquark distribution to be smaller than the corresponding 
longitudinally polarized one, since the latter is parametrically suppressed in the soliton 
size. This expectation is indeed borne out by the numerical results, as will be shown in 
Section [O. 

The behavior of the isosinglet polarized distribution in the limit of large soliton size 
can also be studied using gradient expansion. The gradient expansion is readily derived 
from the spectral representation of the distribution functions, Eqs. ( 13.43]) and Eq. (|3.44| ), 
by formally expanding the quark propagators in powers of derivatives of the pion field, as 
discussed in Section |3~4| , cf. Eq. (|3.50| ). We shall not quote the lengthy expressions here. 
Rather, we shall discuss the properties of these distributions at the hands of the exact 



numerical results in Section 0~2 



4.2 Numerical results 

We now turn to the numerical evaluation of the transversity distributions in the chiral 
quark-soliton model. The numerical calculations are based on the expressions of the 
transversity quark- and antiquark distributions as sums over quark levels in the classi- 
cal pion field of the soliton, Eq. (|3.34j) for the isovector, and Eq. (|3.46| ) for the isosinglet 
distribution. In fact, in the actual calculations we prefer to use the representations of 
the distributions as sums over either occupied or non-occupied levels, Eqs.( p.37|) and 
( |3.38|) viz. Eqs. (|3.48|) and ( p.49|) ; the check of equivalence of summation over occupied or 



non-occupied states offers a very powerful check of the numerical procedure. 

Numerical method. The numerical method for evaluating the sums over single-particle 
quark levels has been described in Ref . ; see also Ref . . The spectrum of quark levels 
is made quasi-discrete by placing the soliton in a spherical 3-dimensional box, imposing 
so-called Kahana-Ripka boundary conditions []55| . The eigenvalues and eigenfunctions of 
the Dirac Hamiltonian in the background pion field are then found by numerical diago- 
nalization, using the eigenstates of the free Hamiltonian (zero pion field) as a basis. One 
can then compute the sum over quark levels, using the fact that the matrix elements of 
the single-particle operators appearing in Eqs.( |3.34"D and ( |3.46| ) between the basis states 
are either known or can easily be computed numerically (see Ref. ||25|| for details). 

Since the Hamiltonian is invariant under combined spin-isospin rotations, and the 
basis states are eigenstates of the sum of the total angular momentum and isospin, it 
is advantageous to convert the expressions for the distribution functions, Eqs. (|3.34|) and 



Q3.46p , to a spherically symmetric form before performing the sums over levels. In this way 
one greatly reduces the number of non-zero matrix elements of the relevant single-particle 
operators between basis states. This symmetrization is achieved by replacing in Eqs. (|3.34j) 
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and ( |3.46 ) the 1- and 3- vector components of the momentum operator, gamma matrices, 
and isospin matrices, by components along two orthogonal 3-dimensional unit vectors, n 
and m, and averaging over their orientations, with the constraint n • m = 0. The matrix 
elements of the resulting "spherically symmetric" single-particle operators between basis 
states of angular momentum plus isospin can then easily be computed. 

In the calculation of quark- and antiquark distributions in the chiral quark-soliton 
model an additional complication arises due to the fact that the relevant single-particle 
operators are discontinuous functions of the single-particle momentum and energy opera- 
tors, which leads to problems when performing the sums over of quark levels in a discrete 
basis. It was shown in Ref. |25| that this problem can easily be circumvented by applying 
Gaussian smearing in the variable x to expressions for the distribution functions. We 
shall apply this method here, using a smearing width of 7 = 0.1, see Ref. for details. 

In the numerical calculations we have used the soliton profile determined in a self- 
consistent minimization of the soliton energy calculated with a Pauli-Villars cutoff 



In the case of the isosinglet unpolarized distribution this choice of profile, combined with 
a corresponding Pauli-Villars regularization of the distribution functions, allowed to pre- 
serve the momentum sum rule for the flavor-singlet quark plus antiquark distributions. 
In the case of the transversity distributions, as shown in Subsection p.4| , the expressions 
for the distribution functions are ultraviolet finite, and there are no sum rules linking 
them to any quantity requiring regularization. For reasons of consistency in the choice 
of model parameters we compute also these distributions with the soliton profile, nucleon 



mass, and moment of inertia, calculated with the Pauli-Villars regularization of Ref. [56 



Although the transversity distributions are ultraviolet finite, numerical calculations 
require that we first evaluate the sums over levels, Eqs. (|3.34|) and (|3.46| ), applying some 
smooth cutoff for levels with large energies, which at the end of the calculation is removed 
by extrapolation to infinity. When performing the numerical calculations in a finite-size 
box, it is important to take into account the "anomaly-type" phenomenon, described in 
detail in Appendix |3.5| . The nature of this phenomenon is non-commutativity between 
the limit of infinite energy cutoff and the chiral limit. In the strict chiral limit, where 
the soliton profile falls off like P(r) ~ 1/r 2 for r — > 00, the sums over occupied and 
non-occupied quark levels would no longer be equivalent. For this reason it is essential 
to perform the numerical calculations in the box with a soliton profile of finite range; for 
example, with a profile falling off as exp(— M n r)/r for r — > 00, where M n is a parameter 
(not necessarily equal to the physical pion mass), and extrapolate to M n — > at the 
end. The condition which must be satisfied for a proper calculation in a finite box is 
M W D ^> 1, where D is the box radius (see Appendix |3.5|) . If this condition is satisfied 
one observes equivalence of the results of summation over occupied and non-occupied 
states in the numerical calculations. In the actual calculations we have used a box size of 
D = (20... SC^M" 1 and values of M w in the range (0.3 . . . 0.6)M. 

Numerical results. The numerical results for the isovector transversity distributions 
are shown in the upper plot of Fig.|l|. It is interesting to first compare the relative con- 
tributions of the bound-state level (dashed line), and of the Dirac continuum (dotted 
line) of quarks, to the quark- and antiquark distributions. Fig.|l] shows that the bound- 
state level gives the numerically dominant contribution to both the quark and antiquark 
transversity distributions. This is in agreement with the observation that the isovector 
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distribution is suppressed in gradient expansion, since the gradient expansion expression 
may be regarded as an estimate of the Dirac continuum contribution. Our results jus- 
tify the approximation made by some of us in Ref. |Q, where only the level contribution 
to the isovector transversity distributions was retained. FigfJ (right column) shows the 
isovector quark- and antiquark distributions multiplied by x. The antiquark distributions 
are shown separately in Fig. |3|. Note that our approach predicts a definite sign for the 
isovector antiquark distribution, 5u(x) — 5d(x) < 0. In the numerical calculation at the 
physical soliton radius this happens because of the dominance of the level contribution, 
which has positive sign, see Fig.[l| The gradient expansion, Eq. (|3.51|) , which becomes 
exact in the limit of large soliton size, also predicts a negative sign of 5u(x) — d~d(x). As 
it is insensitive to the soliton size, the sign of the isovector antiquark distribution can be 
taken robust prediction of this model. 

The different contributions to the isoscalar transversity distribution are shown in the 
lower plot of Fig.[TJ. One sees that in the isoscalar case the dominance of the level contri- 
bution is even more pronounced than in the isovector case. In particular, in the isoscalar 
antiquark distribution the level and continuum contributions all but cancel, leaving the 
total isoscalar antiquark distribution to be much smaller than the quark distribution (by 
a factor of the order of l-2xl0~ 2 , the precise value depending on x). 

Comparing with the longitudinally polarized distributions. It is interesting to compare 
the numerical results for the transversity distributions with those of the longitudinally 
polarized distributions |39|, which are shown in the right column of Fig. 0. One 
notices that, generally speaking, the quark distributions are of similar shape, although 
of different magnitude. The ratios of the numerical results for the transversity to the 
longitudinally polarized quark distributions are well described by the forms 



Su(x) — Sd(x) 
Au(x) - Ad(x) 

5u(x) + 5d(x) 
Au(x) + Ad(x) 



1.25, (4.6) 
2.0 -1.5a;. (4.7) 



These fits apply for values of x > 0.1. (Note that our model, at the present level of ap- 
proximation, cannot be applied to study the small-x behavior of the parton distributions 
at the low scale. Anyway, it is the intermediate- and large-x region of the distributions 
at the low scale of the model which determines the behavior of the distributions in the 
small-x region at experimentally relevant scales due to perturbative evolution.) Note that 
this parameterization is a purely numerical fit. We find it useful to represent the results in 
this way, since we believe the ratios of distributions to be less model-dependent than the 
distributions themselves. In our estimates of spin asymmetries in Section |5| we shall use 
these ratios, together with a parameterization of the longitudinally polarized distributions 
obtained from fits to inclusive DIS data. 

Comparing the transversity and longitudinally polarized antiquark distributions we 
see that they are very different. This is in accordance with the fact that the two isovec- 



tor distributions appear in different orders of the gradient expansion, see Section |3.4| . It 
would thus not be useful to parameterize their ratio, the more since the standard param- 
eterizations for the longitudinally polarized distributions exhibit large differences in the 
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polarized antiquark distributions^]. When computing observables one should rather use 
the results for these distributions directly. 

Results for tensor charges. The tensor charges can be obtained by integrating the 
numerical results for the transversity distributions, or directly by computing the matrix 
elements of the local chirally-odd operators, Eq.( |2.7|) ; both calculations give identical 
results. Our numerical results for the isovector and isoscalar tensor charges are 

5u-5d = 1.06, (4.8) 
5u + 5d = 0.63. (4.9) 

These numbers refer to the low normalization point of (9(600 MeV); however, the scale 
dependence of these quantities is known to be weak. Our results are consistent with those 
of previous calculations in the chiral quark-soliton model |[47|| , which were based on a 
different type of ultraviolet regularization (proper time regularization) . Furthermore, we 



observe a good agreement with the QCD sum rule calculations of Refs. [f43|, 46], and with 



the lattice estimates of Ref.[44 



Discussion of previous calculations of transversity distributions in the chiral quark- 
soliton model. The transversity distributions have been computed by Gamberg et al. in 



a description of the nucleon as a chiral soliton of the Nambu-Jona Lasinio model |42j| , 
which is largely equivalent to the chiral quark-soliton model, and by Wakamatsu and 
Kubota [27[] in same approach as the one used here. In the calculation of Ref.|42j only 



the contribution of the discrete bound-state level was taken into account. It is known 
that in general this simplification, which is not warranted by any parametric limit of the 
model, leads to a number of inconsistencies (e.g. the positivity of the isoscalar unpolarized 
antiquark distribution is violated), as has been discussed in detail in Refs. P4], E5|. In 
the transversity quark distributions it so happens that the contributions from the Dirac 
continuum of quarks are numerically small, as shown by the results of our calculation in 
Fig. [I]. In the antiquark distributions, however, the continuum contributions are seen to 
be numerically important. 

In Ref. a calculation of the transversity distribution was presented including also 



1/TVc-corrections to the isovector distribution, 5u(x)—5d(x). As was noted also in Ref. ||42j| 



these corrections are afflicted by ordering ambiguities of the collective operators, and it 



is not clear if the ordering adopted in Ref. [27] is correct. Short of a satisfactory answer 



to these questions we have chosen to limit ourselves to the contributions appearing in the 
lowest non- vanishing order of 1/N C . A detailed investigation of the i/iV c -corrections will 
be reported elsewhere. 

An important difference of our results to those of Ref.f27j| concerns the ultraviolet 



regularization. In Ref.]p7|1 the Dirac continuum contributions to all distributions was 
subjected to a Pauli-Villars subtraction, although our analysis shows that these quan- 
tities are ultraviolet finite and do not require regularization. In contrast, we have left 
these distributions unregularized. A more detailed comparison with the results of Ref. [T27|] 
is, unfortunately, not possible, since it is not clear how the anomaly-type phenomenon 
discovered in our calculation (see Section ISTST) affects the results of Ref. IE7H, where the 



5 The flavor asymmetry of the polarized antiquark distributions, Am — Ad, was assumed to be zero in 
most parameterizations, e.g. 
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equivalence of summation over occupied and non-occupied levels in the final results was 
not convincingly established. 



4.3 Inequalities in the large N c — limit 



The formal limit of large N c allows one to derive certain positivity bounds for the transver- 
sity distributions, as well as generalizations of the Soffer inequalities, which are stronger 
than those which hold in QCD at N c = 3 [ 40 1 . It is interesting to check if these inequalities 



are satisfied by the results of the numerical calculations. 



In Ref.[33| the following positivity constraints were derived in the large-iV c limit of 
QCD: 



u(x) + d(x) 
3 

u(x) + d(x) 



> \Au(x) - Ad(x) 

> \Su(x) — 5d(x)\, 



(4.10) 
(4.11) 



and similarly for the corresponding antiquark distributions. In the chiral quark-soliton 
model the second inequality is trivially satisfied in the limit of large ultraviolet cutoff, since 
the isoscalar unpolarized distribution is logarithmically divergent, while the transversity 
distributions are finite. For finite cutoff, however, the inequalities could in principle be 
violated by the ultraviolet regularization. In Fig. f| we have plotted the results for both 
the L.H.S. and the R.H.S. in the chiral quark-soliton model. (The isoscalar unpolarized 
distribution has been taken from Ref.[56|.) One sees that the numerical results respect 
the inequalities, both for the quark and antiquark distributions, to a very good extent. 
The small numerical violation should be attributed to the ultraviolet regularization. It is 
interesting that for the quarks we observe saturation of the inequality for the transversity 
distribution, while the longitudinally polarized distribution is smaller than its bound, 
while for antiquarks the situation is opposite: The longitudinally polarized distribution 
saturates the inequality, and the transversity distribution falls short of its bound. This 
qualitative behavior could be a useful guideline for parameterizing the distributions in 
analysis of experimental data. 

Similar is the situation with the "large N c version" of the Soffer inequalities. In Ref. [f|(| 
the following large-iV c inequalities were proven: 



u(x) + d(x) 



+ (Au(x) - Ad(x)) 



> \Su(x) -Sd(x)\ 



(4.12) 



In Fig. H we plot the L.H.S. and the R.H.S. of this inequality, as obtained in the model 
calculation in the large-iV c limit. The same caveats concerning the ultraviolet regulariza- 
tion apply as in the case of the positivity conditions discussed above. As we can see, the 
model results satisfy the inequality within the expected accuracy. 



27 



5 Transverse spin asymmetries in polarized Drell— 
Yan pair production 

With the numerical results for the transversity quark-and antiquark distributions we can 
proceed to make predictions for observables which would allow to extract these distribu- 
tions from experiment. As explained in the Introduction, the transversity distributions 
cannot be measured in inclusive DIS at leading-twist level. In semi-inclusive DIS they en- 
ter together with chirally odd fragmentation function, which, too, are essentially unknown 
quantities |J. The cleanest way to measure the transversity distribution at leading-twist 
level seems to be Drell- Yan (DY) pair production in scattering of transversely polarized 
pp or pp. We shall therefore concentrate on this process here. 

The cross section for DY pair production is a function of the center-of-mass energy of 
the incoming protons, s— {P1+P2) 2 , and of the invariant mass of the produced lepton pair, 
M 2 = (ki + &2) 2 , which is equal to the virtuality of the exchanged photon and where ki/2 
are the momenta of the detected leptons. At the partonic level this process is described by 
the annihilation of a quark and an antiquark originating from the two protons, carrying, 
respectively, longitudinal momenta X\P\ and 3:2^2 -Q The momentum fractions are given 
by Xi/ 2 = (Q 2 /sY^ 2 e ±y with rapidity y — |ln ^1+^) , We consider the case that the two 
protons are transversely polarized relative to the beam direction. In leading-order QCD 
the transverse spin asymmetry of the DY cross section is given by 

^ fcs , M *) = y%(^)w;> , ( 5,) 

J2 f ej q f {x 1 ,M 2 ) q f -(x 2} M 2 ) 

where the sum runs over all species of light quarks and antiquarks in the two nucleons, 
/ = {u,u,d,d, . . .}, The relevant scale here for the parton distribution functions is the 
virtuality of the photon, M 2 . Note that this asymmetry is sensitive to the antiquark 
distributions. We neglect strange quark contributions here; since they always enter in the 
form of a product of a strange quark with a strange antiquark distribution they can be 
expected to be very small. 

In the case of DY pair production in pp rather than pp collisions, using charge conju- 
gation, 5qf/p(x) = 5qjip(x) and qf/ p (x) = qj/ p (x), the above expression changes to 

£/ e j Qf{ x u M ) 9/(^2, M ) 



The transverse spin asymmetry, Eqs.( |5.lD and (|5.2j ), requires the quark- and an- 
tiquark distributions of the individual quark flavors Su(x),5d(x),5u(x) and 5d(x). In 
our approach, based on the large- N c limit, we have computed the isovector distribution 
[Su(x) — 5d(x), 5u(x) — 5d(x)}, which appears in leading order of the l/iV c -expansion, and 
the isoscalar [5u(x) + 5d(x), Su(x) + Sd(x)}, which appears in next-to-leading order, both 
in the lowest non- vanishing order of 1/N C . From these results one should not, strictly 
speaking, recover the distribution of the individual flavors by adding and subtracting the 

6 For questions concerning the reconstruction of the partonic initial state from the event data, see e.g. 
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isovector and isosinglet combination, since there can be l/iV c -corrections to the isovec- 
tor distribution, of the same order in 1/N C as the isosinglet distribution, which are not 
included. Such l/iV c -corrections were studied in Ref . ]27J , however we did not compute 
them here but rather take a pragmatic stand. We directly reconstruct the individual 
flavor distributions from the results of the model calculations. For the quark transver- 
sity distributions we take the model result for the ratios of transversity to longitudinally 
polarized distributions, Eqs.Q4.6D and (|4.7| ), together with the GRSV 95 LO distribution 
for the longitudinally polarized distributions. For the antiquark transversity distributions 
we neglect the contribution of the isoscalar distribution as it is much smaller than the 
isovector one (see Figs.|l| and 0). 

In order to get an impression of the dependence of the prediction for the transverse spin 
asymmetries on the choice of transversity distributions we compare the above results with 
the asymmetries calculated under the assumption that Sq(x) = Aq(x),5q(x) = Aq(x), 
using again the GRSV 95 parameterization for Aq(x) and Aq(x). (This choice, which is 



consistent with the Soffer inequalities, has frequently been made in the literature M, 62 



The transverse spin asymmetries obtained with the two "scenarios" for the transver- 
sity distributions are shown in Fig. |6|. One sees that the differences in the pp asymmetries 
are quite sizable, in particular in the region of small rapidities. Note, however, that these 
observables depend very sensitively on the small antiquark distributions, which may be 
affected by l/iV c -corrections. Even greater differences are seen in the asymmetries for pp 
reactions. Here the asymmetries are dominated by the products of quark distributions in 
the proton, while the contributions from the antiquark distributions are negligible. The 
differences between the asymmetries calculated with our model distributions and those 
with 5q(x) = Aq(x) essentially reflect the numerical enhancement of the transversity 
quark distributions over the longitudinally polarized ones found in the model calcula- 
tions, cf. Eqs.( |4.6p and (|4.7| ). To summarize, we find that our model distributions result 
in significant deviations from what is obtained with that approximation. Unfortunately, 
recent studies, taking into account the limited detector acceptance, suggest that mea- 
surements at RHIC are unlikely to be able to discriminate between the two scenarios 

m. 



6 Conclusions 

In this paper we have studied the transversity quark-and antiquark distributions in a 
dynamical model of the nucleon based on the large- iV c limit. Let us briefly summarize 
the main qualitative conclusions. 

Comparing the transversity and longitudinally polarized distributions, we found that, 
generally speaking, the quark distributions (both isovector and isoscalar) are comparable 
in magnitude. The corresponding antiquark distributions, on the contrary, are very dif- 
ferent. Appealing to the gradient expansion, which, as we saw, gives a realistic numerical 
description of the antiquark distributions, we were able to explain the smallness of the 
isovector transversity antiquark distribution relative to its longitudinally polarized coun- 
terpart on grounds of the constraints imposed on the low-energy effective dynamics by 
chiral symmetry. A measurement of the transversity antiquark distributions would thus 
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be a sensitive test of the role of chiral symmetry in determining the parton distributions 
of the nucleon at a low scale.f] 

The strong differences between the longitudinally and transversity antiquark distribu- 
tions should be taken into account when making predictions for observables sensitive to 
the antiquark distributions. Our estimates show that these differences have a noticeable 
effect e.g. on the spin asymmetries in polarized Drell-Yan pair production. It remains 
to be seen if these asymmetries can be measured to an accuracy that would allow one to 
discriminate between the different predictions |66f . 

Also, we have verified that the quark/ antiquark distributions obtained in the chiral 
quark-soliton model satisfy the Soffer inequalities, as well as the large-iV c inequalities 
derived in Ref.[[|0]. Since these relations could in principle be violated by the ultraviolet 



cutoff required in the model calculation, their fulfillment should be seen as another piece 
of evidence in favor of the Pauli-Villars regularization scheme adopted here [fj4], [25|, |56| 



We have noted a curious anomaly-type phenomenon in the tensor charge and the 
transversity distribution functions. Here we have described this phenomenon using the 
language of sums over quark single-particle in the background pion field specific to our 
large-iV c picture of the nucleon. It is possible, however, that the anomaly described here 
has meaning also outside of the context of this model. To clarify this one should see if 
the observations made here could be stated in a more general, field-theoretic language. 
The analogy with the Fujikawa formulation of the U(l) anomaly, whose field-theoretical 
description is well-known, should be a useful guideline. 
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A Anomaly— type phenomenon in the transversity dis- 
tribution function 



In the calculation of the isovector transversity distribution function in the effective chiral 
theory we encounter an anomaly-type phenomenon: The distributions obtained by sum- 
ming the contributions of occupied and non-occupied quark levels in the soliton are differ- 
ent if the soliton field does not fall off faster than 1/r 2 at large distances. In Section |3l| we 
computed the anomalous difference in the first moment of the isovector distribution, i.e., 
the isovector tensor charge. Using a somewhat different approach one can easily compute 
the anomalous difference in the ^-dependent distribution function. Since the anomaly 
phenomenon in the transversity distribution is of principal theoretical interest, as well as 
of practical importance for the numerical calculations (see Section [4. 2}), we derive here the 



expression for the anomalous difference of the isovector distribution function. 

The representation of the isovector transversity distribution as a sum over occupied 
quark levels, Eq. (j3.37 ), can be written as an integral over a continuous energy variable, 



uj, in the form 8 



[Su(x) - Sd(x)] = / dwp(u) (A.i; 



where the integrand p(u) is defined by the functional trace of the quark "density of states" . 
5(uj — H), with the relevant single-particle operator 



2N c M N q 



5{oj - H)5(cu + P 3 - xM N ) 1+1 ° 1 % 5 -f 1 r l 



(A.2) 



Here, P l denotes the single-particle three-momentum operator, and H the Dirac Hamilto- 
nian, Eq. ( |3.5| ). Similarly, the representation of the distribution as a sum over non-occupied 
quark levels, Eq.( p.38[) , can be written in the form 



[8u{x) - M(z)] non _ occup . = J dup(uj). (A.3) 
The difference between the two representations is thus given as an integral over all energies: 



[Su(x) - Sd{x)] occw _ - [5u{x) - 5d{x)] non _ occu „ = [ dup{u). (A.4) 



We have introduced here an explicit energy cutoff, uq, in order to specify the limiting 
procedure leading to the anomaly. In order to analyze the difference, Eq. ( |A.4| ), in the 



limit of large u>o we write the delta function of uj — H as the imaginary part of the quark 



propagator |25|j: 



[5u(x) - Sd(x)] occup - [5u(x) - $d(x)} non _ oant> 



8 Throughout this section it will be understood that 5u{x) etc. denotes the quark distribution for x > 0, 
and minus the antiquark distribution at x < 0, cf. Eq.(pT 
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ANM 



N 



Im 



du 

— Sp 



H-u-iO 



5(u + P 3 



m ^1±tV 1 1 

xM N ) 757 r 



(A.5) 



We now proceed as in the case of the tensor charge in Section |3~5| , and expand the integrand 
for large uj. This can be done by writing the quark propagator in the form 

1 = H + " (A . 6) 

H — uj — iO H 2 — uj 2 — iOsigncj ' 

and substituting the formal series expansion of this expression in powers of derivatives of 
the pion field, Eq. (|3.50|) , collecting all terms contributing in a certain order of 1/uj. In 
our case the leading contribution in 1/w comes from term with n = 1 in Eq.( |3.50 ); the 
term with n = cannot produce a 7 1 Dirac matrix needed to compensate the 7 1 matrix 
in the operator in Eq. ( |A.5| ) . Taking into account also 75-parity one finds that the leading 
contribution at large uj is given by 

[5u(x) - Sd{x)) occup - [5u(x) - 5d(x)] non _ occup 



AN r M 



N 



Im 



duo . 
2^ 



Sp 



{uj - iTTdk) 



-OJQ 



I+7V ! ! 

x 7; 757 t 



\—d 2 + M 2 — uj 2 — iOsignu;] 



r5(o; + P 3 - xM_ 



(A.7) 



We can evaluate the functional trace in the basis of momentum eigenstates, Sp [. . .] = 
/ d 3 p/ (27r) 3 (p\ ■ ■ ■ \p), inserting complete sets of intermediate position eigenstates in which 
the pion field is diagonal. In addition, we have to take the trace over Dirac and flavor 
indices. In this way we obtain 

[5u(x) - 5d(x)] occnp - [5u(x) - 5d(x)] non _ occup 



N 



Im 



duo 
2^ 



-CJ0 



d 3 p 
(2nY 



Tr 



Dirac 



(uo + 7°7 fc p fc )5(to' + p 3 - xM N ) 



|p| 2 + M 2 



iOsignu;] 



i + 7 y 1 

x 7; 757 



rf 3 xTr flavor l-tM^dkU^i^r 1 



(A.i 



16N r M N M 



d 3 xTr 



flavor 



{ir% 



xlm 



duo f d 3 p {uj + p 3 )5{uj + p 3 



-^ii(x)]} 

xM N ) 



2tt J (2tt) 3 [| p |2 + M 2 - uj 2 - iOsigna;] 



(A.9) 



It remains to compute the integral over the energy, uo, and the momentum, p, in the last 
expression. Instead of computing the integral first and then taking its imaginary part it 
is convenient to take the imaginary part "under the integral" , substituting 

1 



Im 



(uo + p 6 )5(uj + p 3 - xM N )- 



|p| 2 + M 2 — uj 2 — zOsigncj) 



(uo + p 3 )5{u + p 3 - xM N )7c(-signuj)5'{\p\ 2 + M 2 — uj 2 ) 

xM N 5(uj + p 3 - xM N )Tr(-signoj)5' \-xM N (uo - p 3 ) + |p x | 2 + M 2 } . (A.10) 
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In the last step here we have made use of the condition imposed by the first delta function 
in order to simplify the integrand. The integral over p 3 can be taken using up the first 
delta function: 



duo 
~2n~ 



-U>0 



d 3 p 

doj 
2^ 



Im [. . .] 



-U>0 



d 2 p A 
W) 



-xM 



-signcu)5' -xM N {2uo - xM N ) + \v \ + M . (A.ll 



The integral over uo requires some care. We have 
r duo 



^(-sign^)5' -xM N (2uo - xM N ) + |p^| 2 + M 2 

Z7T L 



(A.12) 



+ / 1—5' \-xM N (2uo - xM N ) + Ip^l 2 + M 2 
2n L 



2xM N 2k 



-5 



xM N (2io -xM N ) + \p ± \ 2 + M 2 + (io -> -u ). (A.13) 



The remaining integral over the transverse component of p can easily be performed after 
replacing / — ► / . It gives rise to a sum of two step functions, which in the limit 

oo, are non-zero for values of x in the range 



(2tt)2 

of large cutoff, oo 



- 2u /M N < x < 2uo /M N . 
Collecting everything we obtain from Eq.( |A.9|) 



(A. 14) 



[Su(x) - 6d(x)] occup - \5u{x) - Sd(x)} 



NrMtfM „ /-2a;, 



67T 2 







< x < 



2loc 



non— occup. 



M N ^~ ~ M N JJ L" —J' ( A ' 15 ) 

Eq. (|A.15 ) represents the result for the "anomalous difference" between the sums over 
occupied and non-occupied levels. 

As in the case of the tensor charge, Eq.( |3.57f ), the anomalous difference of the dis- 
tribution functions, Eqs.( [A.15D , is given by a total derivative of the pion field, which is 
non- zero only if the field drops like 1/r 2 at large distances, r — > oo. In this case the 
unitary matrix at large r takes the form 



U, 



cl 



1 + 



kk-~ 



rl[n T 



(A.16) 



where the constant k is related to the nucleon isovector axial coupling constant in the chiral 
limit by Eq. (|3.56 ). The integral of the total derivative of the pion field then becomes: 9 

3 



d 3 xTr flavor ir 1 d 1 {U cl -Ul l ) 



d xTrfj avor 



J2ir k d k (u cl -ul 



Lfe=l 



9 Instead of converting the integral of the total derivative of the pion field into a surface integral at 
r = oo, one also may directly differentiate the asymptotic form of the pion field, Eq.(A.16), and compute 
the volume integral. In this case one would find that the derivative of the l/r 2 -term in Eq.( A.16| ) gives 
rise to a delta function at r = 0, whose volume integral reproduces Eq.(A.17). 
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4nr 2 
hm — — Tr flavor 



£nV=(C/ cl -[/ c \) 



i 

. k=l 



Air k, 2gA 
"~3~ 9F2 



(A.17) 



The anomalous difference is thus given by 



[5u(x) - 5d(x)] occup - [6u(x) - Sd(x)] non _ occnp 



, -N C M N M 6 < x < — - . (A.18) 

127r 2 F 2 N { Mn Mn ) ^ > 

The support of this function increases with the energy cutoff, ujq. A similar phenomenon 
was observed in the anomalous difference in the isoscalar unpolarized distribution in 
Ref. p5|| . The moments of this function (with respect to the variable x) thus exhibit 
power divergences for co>o — >• oo. In particular, the first moment is linearly divergent :0 

oo 

dx {[Su(x) - M*)] occup . - [**(*) - M^U-occup.} = N fjT - (A. 19) 

— oo 

This is in qualitative agreement with the behavior of the tensor charge in the limit of large 
energy cutoff derived in Section |3.5| . Note that the coefficient of the linear divergence 
depends on the details of the ultraviolet cutoff used, so one should not expect to find the 
same coefficients in Eq. ( |A. 19]) and Eq. (|3.57|) . 



10 In the large- N c limit the support of the parton distributions is not limited to the — 1 < x < 1, so 
we must integrate over the entire real axis in order to recover the tensor charge as defined by the matrix 
element of the local operator, Eq.(2.7). 
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[bu - bd] (x) 





Figure 1: The isovector (top) and isoscalar (bottom) transversity quark- and antiquark 
distributions obtained from the chiral quark-soliton model. The functions shown here 
represent [Su =F 5d] (x) at x > and — [5u^f5d](— x) at x < 0. Dashed lines: Contributions 
of the bound-state level. Dotted lines: Contributions of the Dirac continuum. Solid lines: 
Total results (sums of bound-state level and Dirac continuum) . 
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x [8m - bd] (x) x [Am - Ad] (x) 




0.5 x 1 0.5 x 1 



Figure 2: The total isovector (top row) and isoscalar (bottom row) transversity and lon- 
gitudinally polarized quark- and antiquark distributions, multiplied by x. Shown are the 
total results (sum of level and continuum contributions), corresponding to the solid lines 
in Fig. \^.) Solid lines: Quark distributions. Dashed lines: Antiquark distributions. 
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Figure 3: The transversity and longitudinally polarized antiquark distributions, see Fig. [|. 
Left: Isovector distributions, x[5u — 5d](x) and x[Au — Ad](x). Right: Isoscalar distribu- 
tions, x[5u + 8d\(x) and x[Au + Ad](x). Solid lines: Transversity antiquark distributions 
Dashed lines: longitudinally polarized antiquark distributions. 
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Figure 4: The large-N c improved positivity bounds for quark distribution functions, 
Eqs. (\4.1Q) and Eqs.( \4- for the quark (left) and antiquark distributions (right). 
Solid lines: x [u + d](x)/3. Dashed lines: \x [5u — 5d](x)\. Dotted lines: x [Au — Ad](x). 
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quarks 



antiquarks 



0.5 - 





Figure 5: The large-N c improved Soffer bound for the quark (left) and antiquark (right) 
distributions. Dashed lines: \x [5u — 5d](x)\. Solid lines: The large N c Soffer bound (x[u + 
d](x)/3 +x [Au — Ad](x))/2. The small violation is due to the ultraviolet regularization; 
see the discussion in the main text. 
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Figure 6: The transverse spin asymmetry, Att, Eq.( \5.1\ ), in collisions of transversely 
polarized protons, in two different kinematical regions: s = (40 GeV) 2 , M 2 = (5 GeV) 2 
(left), and s = (500 GeV) 2 , M 2 = (25 GeV) 2 (right). Solid lines: Asymmetries calcu- 
lated with the quark- and antiquark distributions computed in the chiral quark-soliton 
model, cf. Fig.^. For the quark distributions we used the calculated ratios of transversity 
to longitudinally polarized distributions, Eqs.( \4-b\ ) and (\4-l\) , together with the GRSV95 
parameterizations for [Am — Ad](x) and [Au + Ad](x). Dashed lines: Asymmetries 
obtained assuming that 6q(x) = Aq(x) and 5q(x) = Aq(x) (q = u,d), using the GRSV95 
parameterizations [||] for Aq(x) and Aq(x). 
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Figure 7: The transverse spin asymmetry, A^ T , Eg. ^5.2j , in collisions of transversely 
polarized protons and antiprotons. The solid and dashed lines correspond to the cases 
described in Fig.^. 
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